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Small separation between reactants, not exceeding 10 -8 — 1CP 7 cto, is the necessary condi- 
tion for various chemical reactions. It is shown that random advection and stretching by 
turbulence leads to formation of scalar-enriched sheets of strongly fluctuating thickness 
rj c . The molecular-level mixing is achieved by diffusion across these sheets (interfaces) 
separating the reactants. Since diffusion time scale is oc 77^, the knowledge of probabil- 
ity density Q(j] c , Re) is crucial for evaluation of chemical reaction rates. In this paper we 
derive the probability density Q(r) c , Re, Sc) and predict a transition in the reaction rate 
behavior from 1Z oc \f~Re (Re < 10 4 ) to the high-Re asymptotics 1Z oc -Re . The theory 
leads to an approximate universality of transitional Reynolds number Ret r ~ 10 4 . It is 
also shown that if chemical reaction involves short-lived reactants, very strong anomalous 
fluctuations of the length-scale r\ c may lead to non-negligibly small reaction rates. 



1. Introduction 

Efficiency of chemical reactions and combustion crucially depends upon number of reac- 
tant moles mixed on a molecular scale. Indeed, for a non negligibly small reaction rate, 
the separation between reacting species must not exceed r ~ 10~ 8 — 10~ 7 cto, which is 
the main reason for immense importance of mixing process [1] . Slow diffusion in laminar 
media leads to an extremely poor mixing and, being the most common mixing acceler- 
ator, hydrodynamic turbulence plays vital role in natural and man-made processes like 
heat transfer, chemical transformations, combustion, meteorology and astrophysics. The 
mixing process in turbulent flows involves three main steps: 1. entrainment, creating 
pockets of material B in a turbulent flow enriched by a substance A; 2. advection and 
stretching leading to formation of thin convoluted sheets of the thickness r\ c separating 
the reactants; This process is often related to as 'mixing by random stirring'. 3. molecular 
diffusion across these 'dissipation sheets ' on a time scale n^/D. If chemical reaction 
and turbulent mixing processes are fast - molecular diffusion is the reaction rate - 
determining step ( Dimotakis (2005), (1993)). Recently, investigation of the role of scalar 
dissipation and dissipation sheets has become an extremely active field (see, for example, 
excellent reviews by Bilger (2004), Sreenivasan (2004), Peters (2000), Dimotakis (2005) 
and recent papers by Celani et al (2005), Villermaux et al (2003) and Bush et al (1998). 
) 

1.1. Kolmogorov - Batchelor phenomenology. 
Based on a classic Kolomogorov's cascade concept, we illustrate the main qualitative 
features of the mixing of two liquids (chemical components) A and B. (A quantitative 
dynamic description will be developed below). Consider a turbulent flow of a fluid A. 
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The integral scale of turbulence, separating the energy and inertial ranges, is L and 
the large-scale Reynolds number is Re — u rms Ljv. At a time t = 0, a blob of the 
same fluid enriched with chemical component B is placed in the flow. For simplicity we 
assume the linear dimension of a blob r = L. At this point, the area of the interface, 
separating A and B substances, is 0(L 2 ). The life- (turn-over) time of this eddy is 
To = T e ddy ~ L/u rms Ls £3 , where u rms is the rms value of turbulent velocity. The 
magnitude of the energy flux across scales is £ = a£U^, ms /L with as w 0.8 — 0.9. In 
accord with Kolmogorov's phenomenology, after w To seconds, by nonlinear interactions, 
this eddy is transformed into another one of the linear dimensions L x% \ = L/2 and 
Ly.x = L z ,i rts \[2L. Since the gradient in the x-direction is largest, the characteristic 

2 2 

time-scale of this, "daughter" structure is: T\ w ^-f-2 _ 3 . Then, after n»l steps 

2n hi I? 

L x ,n — 2 n L; T n ~ 2 3 — 7- ~ 2 " — (1-1) 

£3 v 

where the viscous time is denoted as r„. We see that the viscous time of a structure 
strongly decreases with the number of 'cascade steps'. The time needed to form the 
smallest Kolmogorov eddies on the scale t\k ~ Ti?e~* is thus: 

T K = T + Ti + • • -T„ « 2.7T « §i (1.2) 

£ 3 

where n « | " ^ e w lni?e ^> 1. It is important that the Reynolds number Re = u rms L/v 
used in the above relations is to be distinguished from Reu = UL/v where U is the mean 
velocity in the flow. Typically, numerically, Re < Reu. The fact that Kolmogorov's scale 
r\K and dissipation rate £ are formed on the time-scale of a single large-scale turn-over 
time (T e ddy) has been tested in various numerical experiments. 

If the Schmidt number Sc = u/D w 1, in accord with the cascade picture, the Kolmogorov 
scale is the smallest legth-scale created by turbulence and at times t > Tk, the mixing 
proceeds by molecular diffusion. 

If however, Sc 3> 1 and the scalar diffusion is extremely inefficient, then, after formation 
of the Kolmogorov scale tjk, the stretching process by the large-scale (r << t]k) velocity 
field leads to generation of the ever thinner scalar-enriched sheets until f] c (t) ~ t\b ~ 
t]k/VS~c <C r\K- It is only after that, the scalar diffusion takes over the mixing process. 
The mean width of these sheets is: t\b ~ tIk/VSc, where j]b is called the Batchelor 
(1959) scale. At this stage, the substances A and B, separated by the distance t/b, can 
mix only by molecular diffusion on a time scale Td ~ V%/D ~ t \kI v ~ —Re~^ ■ 
To get a feel for "real life charm", we consider a simple numerical example. In gases 
where Sc » 1, the length-scales tjk w t\b- In liquids, the situation is different and 
Sc w 600 - 3000 (Dimotakis (2005), (1993)) which means that t] K /t] B w 25 - 50. 
Thus, while the scales t)b and r\K can be very different, the corresponding diffusion and 
viscous times scales Td ~ Vb/D ~ t \kI v ~ T v are of the same order. Moreover, it has been 
shown in various numerical and experimental studies (see for example Monin and Yaglom 
(1975)) that the relevant viscous scale 77^ is numerically larger than Kolmogorov's one: 
T} v = a v rjK with a v 10 — 20. Thus, the diffusion time is, in fact, Td = rjl/v w r„. (The 
origin of the large factor a v will be discussed below.) 

1.2. Low and large Reynolds number behavior of reaction rates. 

If Sc » 1, the process consists of two steps: 1. Formation of structures (sheets) on 
the viscous scale r\ v w a„LRe~^; 2. Further stretching of the scalar fieled toward n c ~ 
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Ci v r\K I '\/Sc. It can be shown readily (see Monin and Yaglom (1975)) that, after initial 
formation of the dissipation scale r\ v , the distance between two particles across the sheets, 
stretched by the large-scale velocity field decreases with time as r{t) = r\ v exp(— -ft). Thus, 
the scalar dissipation scale is formed on a time - scale r c given by the relation: 

r\ v e~ 1T ' w n^/VSc 
where 7 = a 7 ^/f with > a 7 > This gives: 

t c = - — — Re~^ In Sc 

The diffusion time across scalar dissipation sheets is: 

d*,L 2 Re~§ , L „ 1 
T d w — w < fle 3 

n rrns 

Comparing the above relations we come to a non-trivial conclusion: although the charac- 
teristics times scale with the Reynolds number as r c oc Td oc \/Re, the stretching process 
leading to both sheet thinning and increase of the scalar gradient competes with the si- 
multaneous concentration-gradient-decreasing diffusion across the developing interfaces. 
However, due to the large magnitude of numerical factor a„ w 10 — 20, the process of 
the sheet-thinning by the large-scale stretching is numerically much faster than diffusion 
across the sheets. This means that during the stretching stage, the scalar diffusion across 
the interfaces can be neglected. The factor a v w 10 — 20 will be calculated below and it will 
become clear that its large magnitude is a consequence of complex small-scale dynamics 
of intermittent turbulence. 

Thus, the ratio of the mixing (inviscid) and diffusion times can be estimated as: 

Teddy 



Re 



Td af, 

and T e ddy /rd ~ 1 for Re w a 2 v w 10 4 . 

In the flows with Re < 10 4 , molecular diffusion across the dissipation sheets is the 
longest, rate-determining, process. In the wall flows, this Reynolds number corresponds 
to Re v < 2 - 3 x 10 5 . It is only when Re v > 10 5 - 10 6 , the rapid diffusion through 
extremely thin interfaces, is dynamically irrelevant for the mixing process and inviscid 
mixing time T e ddy is the rate-determining step. Based on the above considerations, one 
can expect a transition from the diffusion- to advection - dominated mixing at 

Re tr w a 4 (1.3) 
If mixing is the reaction rate ( 1Z ) determining process, we expect 



V V 3 / 

U OC -5- W -r^Re^ = TeddyVRe Re < Re tr 

V 2 k L 2 

n oc w ^Re 1 = TeddyRe Re > Re tr (1.4) 

In accord with (1.4), the transitional Reynolds number depends only upon coefficient a u , 
characterizing small-scale properties of turbulence. Thus, we can conclude that, since the 
small-scale property of turbulence a v s=s 10 — 20 is a more or less universal number, inde- 
pendent upon type of the flow, the derived Re tr ~ 10 4 must be approximately universal. 
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A mixing transition leading to the Reynolds - number - independent reaction rate at 
approximately universal Reynolds number Rctr ~ 10 4 , has been observed in experiments 
by Dimotakis (2005), (1993). 

2. Statistical Description of Disipation Structures. Evaluation of a v . 

The Kolmogorov theory (K41), treating the u.v. cut-off tjk ~ LRe~i = const, completely 
disregarded the non-trivial dynamics of the dissipation range fluctuations. It became clear 
recently that the dissipation scale is not a constant number but a random field defined 
as : 



8 v u \(u(x + rj) — u(x))\ 
where S y Ui = Ui(x + y) — Ui(x). In this form the relation for n was derived (Yakhot 
(2003, (2006)), Yakhot and Sreenivasan (2004), (2005)) from the dissipation anomaly, first 
introduced by Polyakov (1995) for the case of Burgers turbulence and later generalized to 
the Navier-Stokes turbulence by Duchon and Robert (2000), Eyink (2003). Even earlier, 
the fluctuating dissipation scale was used by Paladin and Vulpiani (1987) in the context 
of multifractal theory The local value of the Reynolds number Re ~ u^rj/v w 1, with 
5 v u = u n as a typical speed of an eddy of linear dimension 77, was mentioned in Landau 
and Lifshitz as a criterion for the onset of viscous dissipation, as early as 1959. 
The physical meaning of the dissipation scale 77 is understood as follows. As r — > 0, the 
velocity field is analytic, so that S r u w d x u(x)r. Thus, in this limit S n = (S r u) n oc r n . On 
the other hand, in the inertial range S3 oc r which, contradicting small-scale analyticity of 
velocity field, cannot be valid in the limit r — > 0. Moreover, in isotropic and homogeneous 
turbulence, the dissipation rate E oc v(d x u x ) 2 = — i^lim r ^o -§piS2(r) oc i^lim r ^o r ~ 3 = 
0(1) is regularized by viscosity v — > 0. This leads to the definition of the dissipation scale 

1 — - 3 
t\k '■ £3 oc vr/ K 3 and, since £ w u^, ms /L, to Kolmogorov's estimate for tjk oc -Re". 

Based on the above considerations, defined by (2.1), n is a scale separating analytic 
(r <C rj) and singular (r 3> rj) contributions to turbulent velocity field. The relation (2.1) 
is an order of magnitude estimate of the "dissipation scale" and, in general, r\ ss j^- where 
a is a velocity- field -independent factor which was investigated in numerical simulations 
by J. Schumacher et. al. (2007) . Interested in qualitative aspects of the mixing process, 
we, for now, neglect this coefficient and use expression (2.1). 

According to the theory (Yakhot (2003)), there exist an infinite number of "dissipation 
scales " r\ n separating smooth S n oc r n (r <C r] n ) and singular S n oc r^ n (r >• r) n ) 
intervals of the moments S n — (6 r u) n . This fact has been decisively demonstrated in 
numerical experiments of Schumacher et al (2007). According to the analytic theory 
(Yakhot (2003)), in general 

LRe- f3 <n n < LRe~i 

where (3 ~ 1 and the Kolmogorov scale t\k ~ ?72 is only one of the possible dissipation 

scales. This property of turbulence will be important in what follows. 

Let us consider an equation for concentration c of a passive scalar advected by velocity 

u: 



— + u • Vc = L>V 2 c; V-u = (2.2) 

We assume that the velocity field u is governed by the Navier-Stokes equations and 
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u = U + u' with U and u' corresponding to quasy-regular (sometimes time-dependent) 
and chaotic (turbulent) contributions, respectively. The fluctuations of the scalar dissi- 
pation field N — Z)(Vc) 2 , governed by (2.2), have been investigated in great detail for 
the case of the <5-correlated in time, large-scale, velocity field (Kraichnan-Batchelor prob- 
lem), where the stretched exponential tail of the distribution function P(N) was derived 
in the range N » N ( Chertkov et. al (1998), Gamba et al (1999)). In this paper we are 
interested in statistical properties of "molecular diffusion sheets" , which are not directly 
related to the tails of the scalar dissipation rate distribution. 

The scalar field is analytic, so that for r — > 0, c(x + r) — c(x) ss d Q^ r. In addition, in 
the scalar "inertial range" r ^> ijk (Monin and Yaglom (1975): 



S%' c = {u(x + r) - u{x))(c(x + r) - c(x)) 2 = --Nr (2.3) 



4^ 



where N — D(-j^-) 2 = 0{1). In the case Sc >> 1, there exist an additional scalar 
"rough" range r]B < r < t]k, where Ss(r) oc lnr. It is only at the scales r <C i]b, 
the scalar field is smooth. By definition, the length scale r w n c is the scale separating 
analytic and singular contributions to the scalar field c(r, t). As follows from (2.3), in the 
inertial range, the scalar field c(x) is not differentiable and in the limit D — > 0, one has 
to be careful with evaluation of spatial derivatives of the scalar field. 
From the equaion (2.2) we have 

8c 2 

— + u • Vc 2 = 2DcV 2 c (2.4) 
at 

and introducing the "point - splitting" c(±) = c(x ± y) and u(±) = u(x ± y), derive: 
d(c(+)c(-) 



dt 



(V+u(+) + V_u(-))c(+)c(-) = D{V\ + VI )c(+)c(-) (2.5) 



Equation (2.5) involves derivatives of singular ( in the inertial range) functions. Thus, in 
the limit y — > ; the exact equation (2.4) for the scalar variance can appear from (2.5) only 
if singular and regular contributions balance separately. Indeed, in the limit y — > n c — > 0, 
taking into account that J| = = —g§— and repeating all steps presented in detail 
in (Yakhot (2006)), we derive: 



-^-(SyUiiSyC) 2 ) + 2V+u(+)c 2 (-) + 2V_u(-)c 2 (+) = -ADS y c^S y c (2.6) 

The relation (2.6) is exact locally in space and time. Since in incompressible isotropic 
turbulence the velocity- scalar correlation function Ui(x)c 2 (x') = (Monin and Yaglom 
(1975)), it is clear that averaging over a "ball" of radius rj c (Duchon (2000), Eyink (2003), 
the second and third terms in the left side of (2.6) disappear giving the locally valid 
estimate for the scalar dissipation scale, independent upon specific model of turbulence: 



As follows from this relation, the random variable n c depends upon local values of velocity 
fluctuations. In the most interesting and important case Sc = v/D 3> 1, on the scale 
n c <C f] the velocity field is analytic, giving: 
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Figure 1. Inertial range: r > 770. Probability density P(^r,r) vs 8 r u/r a . a. r=0.01; 

b r=0.01; 0.1; 0.5 



~ "sTT ~ JT~ ~ ~ = (2.8) 

lc du(x) X u y I I V I 

dx n 

Thus, the probability of the scalar dissipation scale Q{j] c ) is evaluated readily from the 
PDF Q(rj) calculated in Yakhot (2006). This result leads to some important consequences: 
If 



Td ~ 1]l I ' v = a lVK I ' v ~ C? v T e ddy / ' VRe (2.9) 

to evaluate the reaction rate and the proportionality coefficient a„, we need the proba- 
bility density function Q(rj c ) = \* r ScQ(rj^fSc). 



3. Probability densities. 

In what follows we set L = 1, so that £ = r < 1 and if the moments of velocity increments 

S n ,o = (S r u) n = A(n)r^ n , then the probability density function can be found from the 
Mellin transform: 

1 />«oo 

P{S r u,r) = — A(n)r^ n \S r u)- n dn (3.1) 

where we set the integral scale L and the dissipation rate £ equal to unity. Indeed, 
multiplying (3.1) by (S r u) k and evaluating a simple integral, gives Sk,o — A(k)r^ k . Under 
different name this transformation has been used in Tcheou et.al (1998) in the context 
of multifractal theory of turbulence . With the Gaussian large-scale boundary condition 
for the probability density at r = 1, the amplitudes A(n) = (2n— 1)!! and, for the values 
of n < 1/6, we can use the Taylor expansion of the function giving w (an — bn 2 ). 
The detailed theory for this case has been recently developed in Yakhot (2006) with the 
result (SrfU = u): 

2 2 (In— ^) 2 

P(u,r) = / e~ x exp[ A r *f\ ]dx (3.2) 

V ' iru^4b\lnr\ J-oo 46|lnr| J V ; 

The PDF P(u/r a 1 r) is plotted on Fig. 1 for a few values of the displacement r. Here 
a rts 0.38 and b = 0.017. On Fig. lb we see the broad tails which, as a sign of strong 
intermittency, cannot be collapsed on a single curve. 

According to its definition, the dissipation scale is a linear dimension of a structure 
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Figure 2. Logarithm of probability densities of the scalar dissipation scale Q(^,Re, Sc) ; 
Re = 10 4 , Sc = 1 (right curve), Sc = 25 (left curve). rj = LRe'T^ with a = 0.38. 

defined by the local value of the Reynolds number Re v = r/S^u/v = O(l). Experimentally, 
probability density Q(t], Re, Re v ) is found by fixing the displacement r = rj and counting 
the events with rjd^u/v = Re v = a. This algorithm has been used by Schumacher (2007) 
in his numerical investigations (see Fig. 3). From the formula (3.1) we have: 

-i /.joe 

P{5 v u) = P{u ri ) = — A(n)v^u-^- n dn (3.3) 

U V J -too 

Fixing L = £ = u rms = 1 gives the large-scale Reynolds number v = 1/Re and taking 
into account that by virtue of (2.1) u n Re/u rms — LS J U ss L/rj gives for the probability 
density Q(jj) (in what follows we denote ^ = rj): 

-I POO POO 

Q{rj,Re) = - e~ x 'dx dne ™in(r,» +1 ^fe)-h„ 2 in, 

V J — OO J —OO 



(3.4) 



and: 



1 f°° _ 2 In^+VggRe) 1 /' 

Q(ri,Re) = — . , , / e x dxe <"> ln " = . / 

vVWtinj i.oo v J4b\ln V \ J- 



1 f°° _„ a In^+Vggge) 1 f°° _ , 1° 2 «^)° +1 ^) 

dxe 4b i '-ii 
(3.5) 



where 770 = LRe^ 1 ^. As expected, the probability density of dissipation scales is ex- 
pressed in terms of the ratio rj = 77/770 and the width of the distribution is the weak 
function of the Reynolds number. The PDFs Q(rj/rio, Re) for Re = 10 4 is shown on Fig. 
2 for two values of the Scmidt number Sc = 1 and c = 25. 



4. Mean dissipation scale and diffusion time. 

Now we can evaluate the moments of the dissipation scale: 

& = {^y - r^YQin, m )dn (4.1) 

Vo Jo Vo 
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and mean diffusion time: 

poo 

- A = rflv=\ V 2 Q(r,, Vo ,Re) dn/v (4.2) 
Jo 

The numerical results slightly vary with position of the maximum of PDF Q (77/770), 
which depends upon the magnitude of parameter a in the expression anu^/v = 1. If 
we choose a so that, in accord with numerical simulation of Schumacher (2007), the 
maximum is set at i] max /i] n w 2, then numerical integration (4.1) yields 77 = a„rj n sa 7?7o 
and T? 2 w I2O770. If however, r] max w 5?7o, we derive 77" = 0^770 ~ 13.5?7o, which is close 
to the outcome of Dimotakis's (2005) physical and Gotoh-Nakano's (2003) numerical 
experiments. In this case, ?7 2 ~ 500r7 O . In general, based on (4.1), (3.5), rj 2 S> a%T} 2 . 
This result leads to important conclusion: Due to strong intermittency of the dissipation 
scales, the 0(rf jv) scalar diffusion time is much longer than w cq,rf K , calculated on 
the basis of Kolmogorov 's phenomenology. Therefore, molecular diffusion, as a reaction 
rate determining process, is even more restrictive than previously thought. We can also 
conclude that these fluctuations are responsible for the "large" magnitude of a constant 
a v w 10-20. 

4.1. Mixing reactants having finite life-time . 

One can define dimensionless Damkollcr number Da = (r c + Td)/r reactio „, which is the 
ratio of hydrodynamic mixing time scale to T reac tion, characterizing the reaction rate 
between perfectly mixed reactants A and B. In case of fast reactions, we are interested 
in here, Da » 1. Below, we will show that in some cases, to describe chemical reactions, 
Kolmogorov's cascade picture, though elegant and illuminating, is not sufficient. 
In what follows we consider a simple example of a model photo-chemical reaction A*+B = 
AB + hv where A* is a component A in an initially prepared electronically excited 
state characterized by a life-time r e . By definition, the finite life-time r e implies time- 
dependence of concentration of excited states ca* = c(t = 0) cxp(— — ). As above, we are 
interested in a diffusion-dominated limit Re < 10 4 . A chemical reaction is possible only if 
diffusion time Td ~ cy 2 r}%/D < r e , a v T)B < V T eD and 77^ < ^/r e v. This gives, in addition 
to the Damkohller mumber, a dimensionless reaction criterion: 



Y = a u ^j^Re-i < 1 (4.3) 

When the life-time r e is very small, according to the relation (4.3), which is the outcome 
of Kolmogorov's phenomenology, the photo-chemical reaction is impossible. Now we will 
show that due to strong fluctuations of the scalar dissipation scale 77c this conclusion 
must be modified. 

Let us consider the advection-diffusion equation for concentration of a passive scalar 
c = ca* undergoing a chemical reaction with another one of concentration cq. If the 
maximum separation, for which reaction is still possible (reaction radius) is A, we can 
define the probability to find a molecule B within a sphere of radius A surrounding the 
reactant A* as ?(r < A) and write the balance equation for the reactant A as: 

^ + u- Vc = DV 2 c- cP(r < A)w AB - — ; V-u = (4.4) 

OX T e 

where wab is the reaction rate of perfectly mixed reactants separated by the distance 
r < A. The 0(l/r e ) term in (4.4) accounts for the finite life -time of one of one the 
reactants A* . It is clear that V(r < A) depends upon concentration cb- As was mentioned 
above, the reaction radius A, depending on the overlap of molecular orbitals, is very small 
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and, during the mixing stage, when r(t) » A, the probability V = and the 0(wab) 
chemical contribution to the balance equation (4.4) can be neglected. The relation (4.4) 
illustrates importance of the molecular-level mixing process in chemical kinetics. 
By substitution ca = c — > ce~~ , the remaining equation is transformed into (2.2). It is 
clear that the reaction rate is not negligibly small only if the mixing time < r e . 

We illustrate the qualitative features of the process on a numerical example. If rj/rjo ~ 
13, (Gotoh/Nakano (2003)), then, based on the PDFs computed above, tH/tq w 500. For 
Te/Td wj<1, the naive (mean- field ) reaction yield, proportional to the concentration 
c oc cxp(— Td/r e ) , is negligibly small. However, defining r e w rjl/u, gives ^ w ^/500g. 
This result means that a finite a fraction 



F= / Q(n, m ,Re)d V (4.5) 
Jo 

of dissipation sheets with n < n e w ^/SOOgTyo does contribute to the nonzero reaction 
rate. Since s=s t„, the integral is evaluated using the probability density of velocity 
dissipation scales. 

Taking, for example, g w 0.02, we see that only the sheets of the thickness r\ < VTOijo w 
3?7o contribute to this reaction. The fraction of the dissipation structures satisfying this 
condition is F m 0.3. If g « 0.01, we find F«5x 10~ 4 . We can see that, due to strong 
fluctuations of the dissipation scale, the reaction is not negligible even when r e <C Tj7. 



5. Conclusions. 

To conclude: the scalar and velocity dissipation scales -n and rj c in turbulent flows are not 
constant numbers but describe random fields with r\ c w -D/ff and rj w v/S^u, respec- 
tively. In an important case Sc 3> 1, these scales are related as r\ c w r]/VSc. Based on 
the Mellin transform and Taylor expansion of the scaling exponents of velocity structure 
functions, the probability density of the scalar dissipation scale has been derived. Two 
main results of this paper are: 1. Due to strong small-scale intermittency, the calculated 
mean thickness of a dissipation sheet is a u rj K where a u m 10 — 13. Extremely strong 
intermittency leads to a long scalar- transport time Td across the sheets and, in the flows 
with Re < aj, rts 10 4 , to diffusion as a reaction rate- determining step . Therefore, the 
reaction rate is: 1Z cx \[Rb for Re < 10 4 , and 1Z oc Re in the interval Re > 10 4 . 2. Even 
when the life - time of the reactants is very short, due to the dissipation scale fluctua- 
tions, the reaction can proceed via diffusion across thinnest dissipation sheets r\ c «rf^- 
In this case, the fluctuations lead to the non-negligibly small reaction rates. This result 
may be of importance for reactions involving short-lived radicals, excimers and other 
cases. We believe that experimental investigation of the sub-Batchelor scale dynamics of 
the mixing process is an extremely interesting and urgent task. 

The theory presented in this paper is based on the dissipation scale definitions (2.1), 
(2. 7), (2. 8), derived from the dissipation anomaly. This algorithm was numerically com- 
pared by Schumacher (2007) with the one based on the isosurfaces of the scalar dissipation 
rate. The obtained PDFs Q{rf), though qualitatively similar, had quite substantial quan- 
titative differences (see Fig. 3). Since the relations based on dissipation anomaly (2.1) 
have been derived directly from equations of motion, we believe they are much better 
justified. 

Interesting and stimulating discussions with N. Peters, U. Frisch, A. Kerstein, J. Schu- 
macher, E. Villcrmaux and K.R. Sreenivasan are gratefully aknowledged. 
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Figure 3. Logarithm of not normalized probability densitiy of the disspation scale 
P(t = rj c ) = Q(j^,Re) for Re ~ 10 3 and a few different values of parameter a. The PDFs 
were numerically evaluated by Schumacher (2007) using two different methods: 1. Algorithm 
based on the dissipation anomaly (eq. 2.1). (upper curve); 2 From the isolevels of dissipation 
rate £. (lower curve). Broken lines denote Kolmogorov's scale. 
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